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ABSTRACT 

We present here a comprehensive study of the optical/near-infrared (IR) upper limits for 
gamma-ray bursts that have an X-ray afterglow. We have extrapolated the X-ray afterglows to 
optical wavelengths based on the physics of the fireball blast wave model (e.g. Rees & Meszaros 
1992; Meszaros & Rees 1997), and compared these results with optical upper limits for a large 
sample of bursts. We find a small set of only three bursts out of a sample of 20 for which the 
upper limits are not compatible with their X-ray afterglow properties within the context of any 
blast wave model. This sparse sample does not allow us to conclusively determine the cause of 
this optical/near- IR deficit. Extinction in the host galaxy is a likely cause, but high redshifts and 
different afterglow mechanisms might also explain the deficit in some cases. We note that the 
three bursts appear to have higher than average gamma-ray peak fluxes. In a magnitude versus 
time diagram the bursts are separated from the majority of bursts with a detected optical/near- 
IR afterglow. However, two GRBs with an optical afterglow (one of which is highly reddened) 
also fall in this region with dark bursts, making it likely that dark bursts are at the faint end 
of the set of optically detected bursts, and therefore the dark bursts likely form a continuum 
with the bursts with a detected optical afterglow. Our work provides a useful diagnostic tool for 
follow-up observations for potentially dark bursts; applied to the events detected with the Swift 
satellite, it will significantly increase our sample of truly dark bursts and shed light upon their 
nature. 

Subject headings: gamma rays: bursts — dust, extinction 



1. Introduction 

After the discovery of the first gamma-ray 
bursts (GRBs) (Klebesadel et al. 1973), it took 
nearly 30 years before emission at other wave- 
lengths (afterglow) was discovered (Costa et al. 
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1997; Van Paradijs et al. 1997; Frail et al. 1997) 
and the distance scale was determined (Metzger 
et al. 1997). We know today that the main rea- 
sons for this delay were the relatively rapid de- 
cay of the afterglows and the inaccuracy of the 
initial burst localizations: (optical) searches for 
afterglows have to be performed within one or 
two days after the burst has been triggered in 
gamma rays and before the afterglow flux level 
has decreased below the telescope detection limit. 
These searches have to be performed with large 
aperture telescopes (requiring arc-minute sized er- 
ror boxes) or with sufficiently long exposure times 
on smaller telescopes. Both options were rarely 
feasible for a long time, as either the delay be- 
tween the occurrence of the GRB and the avail- 
ability of its position was too large, or the size of 
the error box of the GRB was too large to per- 
form efficient follow-up observations. The rapidly 
available and accurate positions provided initially 



by the Wide Field Cameras (WFCs) on board the 
Italian-Dutch BeppoSAX satellite, and later with 
HETE-II, combined with X-ray afterglow counter- 
parts localized with the BeppoSAX Narrow Field 
Instruments (NFIs), Chandra and XMM-Newton, 
made rapid follow-up observations at longer wave- 
lengths possible, with a typical error circle diame- 
ter of the same order as the field-of-view of a large 
optical telescope. 

It was noted very soon, however, that despite 
rapid follow-up and deep imaging, some bursts es- 
caped the detection of an optical counterpart (e.g., 
Groot et al. 1998), and have been traditionally 
designated as 'dark bursts'. Notably, an optical af- 
terglow was found in only 35% of the BeppoSAX 
bursts for which an X-ray (2-10 keV) afterglow 
was detected (see e.g. De Pasquale et al. 2003). 
Interestingly, however, eight out of ten bursts re- 
cently detected with the soft X-ray camera (SXC) 
on board HETE-2, have a detected optical or near- 
IR afterglow, increasing the detection rate to 80%. 
This discrepancy in their dark bursts rates reflects 
rather the current improvement of the commu- 
nity response in GRB alerts, as well as today's 
enhanced observational capabilities, than an in- 
trinsic difference in the properties of the GRBs 
detected by the two satellites. 

Several explanations have been brought for- 
ward to address the absence of detectable opti- 
cal flux, apart from adverse observing conditions: 
dark bursts could be intrinsically faint at opti- 
cal wavelengths (possibly reflecting a difference in 
their central engine from 'normal' GRBs); they 
could be at high redshift, which would result in 
their optical/near-IR light being suppressed due 
to hydrogen (Lyman-a) absorption; or they could 
be heavily extincted by gas and dust. The lat- 
ter would not be unexpected, since the afterglows 
are located close to the center of light of their 
host galaxy (Bloom et al. 2002a). Furthermore, 
it is known now that at least a sub-class of GRBs 
originate from the core collapse of a massive star 
(Hjorth et al. 2003b; Stanek et al. 2003). As 
these massive stars are very short-lived, they are 
expected not to travel far away from their birth 
place, which is most likely a very dusty environ- 
ment. There are cases, however, where the host 
galaxy shows a high column density of neutral hy- 
drogen, as expected for massive-star progenitors, 
but where strong upper limits are set on the red- 



dening (and thus the dust content) in the host 
galaxy; e.g. for GRB 020124 (Hjorth et al. 2003a) 
and for GRB 030323 (Vreeswijk et al. 2004). In 
such cases, it is possible that X-ray dust destruc- 
tion might free the way for the afterglow radiation 
to be still detected (Fruchter et al. 2001). 

It is important to note here that the current 
designation of a burst as 'dark' may reflect either 
that no afterglow has been found by optical/near- 
IR searches, or that the optical counterpart is ex- 
tremely faint initially and it declines very rapidly, 
thus precluding detection when observed later 
than one day after the burst (see e.g. Fynbo et al. 
2001; Hjorth et al. 2002; Berger et al. 2002). In the 
present study, however, we define as 'dark', bursts 
that fall below an upper limit in their afterglow 
flux, which we derive using the fireball model. In 
this physical definition of dark bursts, the input 
for our calculations is obtained from observations 
in other wavelengths (X-ray or radio) where the 
afterglow was detected. Of course, if there is no 
detection at any wavelength, this definition is not 
applicable, and one has to revert to the original 
empirical definition. Thus, to define our sample of 
'physically dark bursts', we have limited ourselves 
to those with a detected X-ray afterglow, from 
which we can predict the optical/near-IR fluxes. 

The typical power-law decline of the X-ray and 
optical afterglows as well as the afterglow broad- 
band spectra obtained from observations at ra- 
dio, millimeter, infrared, optical and X-ray wave- 
lengths, show that the standard fireball model 
proposed as an explanation for these phenomena 
(Meszaros & Rees 1997; Wijers et al. 1997), is to a 
large extent acceptable. This standard model has 
been modified to account for the possible collima- 
tion of the outflow of the gamma-ray burst and 
its afterglow (e.g. Rhoads 1997, 1999), and for the 
expansion of the afterglow into a medium with a 
density gradient (e.g. a stellar wind; Chevalier & 
Li 1999). In general the model can be constrained 
by the available observations of a GRB afterglow. 
On the other hand, if observations in a partic- 
ular wavelength interval are missing, it is possi- 
ble to use the afterglow model to predict the flux 
(as a function of time) in this wavelength regime, 
which is what we have done below. The outline 
of this paper is as follows: in Sect. 2 we describe 
the method used to constrain the optical/near-IR 
magnitudes expected from the X-ray afterglows. 



2 



The results are analyzed in Sect. 3 and discussed 
in the context of the current ideas to explain dark 
bursts. Finally, we show how our results can be 
used as a diagnostic tool for future observations 
and studies of dark bursts. 

2. The method 

2.1. X-ray measurements 

We have collected the X-ray afterglow obser- 
vations from various sources (refereed journals, 
proceedings, GCN Circulars) and obtained a to- 
tal of 20 X-ray afterglows with no optically de- 
tected burst, as well as 1 X-ray afterglow for an 
optically faint burst (GRB 020322). The available 
X-ray fluxes were first converted to the frequency 
corresponding to the observed energy band loga- 
rithmic average, instead of the usually quoted flux 
in the 2-10 or 0.2-10 keV range. This was done 
using e.g. F 2 - W = F av x (^)^dE, where (3 is 
the energy spectral index, E av the logarithmic av- 
erage and F av the corresponding flux (e.g., for the 
2-10 keV interval, E av ~ 5 keV). Since the quoted 
fluxes are generally an average of the flux in the 
observed time interval, we interpolated the flux to 
a logarithmic average t aw using the temporal decay 
index 8 in a similar way. The choice of this spectral 
and temporal average ensures that our fluxes are 
not heavily dependent on the power-law indices (3 
and 8. 

In cases where no spectral index was reported, 
we have assumed a Crab-like spectrum. The er- 
ror induced by this assumption is generally small, 
due to the choice of the logarithmic average as 
our pivot energy: for example, for a F2-10 nux 
of 5 • 10~ 9 crg cm _2 s _1 and (3 between 0.6 to 1.4, 
the F 5 flux varies only between 5.8 • 10~ 10 and 
6.4 • 10~ 10 erg cm _2 s _1 keV~ 1 . We have assumed a 
20% error in the flux and the temporal or spectral 
indices of the X-rays if no error was listed. Table 
1 lists the X-ray afterglows which have been used 
in our analysis. 

2.2. Optical data 

Optical and near-IR data were largely obtained 
from the GCN Circulars, with additions from ref- 
ereed journals and conference proceedings. The 
i?-band upper limit for GRB 000214 was obtained 
during our observing programme on the 2.2m ESO 



telescope at the La Silla observatory. The data 
were reduced in a standard fashion and calibrated 
using a set of Landolt standard stars. 

All optical and near-IR magnitudes have been 
corrected for Galactic extinction using the COBE- 
DIRBE maps (Schlegel et al. 1998). Corrected 
magnitudes have been converted to fluxes using 
the Vega flux values provided by Fukugita et al. 
(1995) for optical magnitudes and those by Beck- 
with et al. (1976) for infrared observations. 

Where the quoted upper limit was the Digitized 
Sky Survey (e.g. "we find no new sources down to 
the limit of the DSS"), we have assumed a limiting 
magnitude of 20. Unfiltered observations were left 
out of the analysis, unless the authors have given a 
corresponding filter and magnitude. Table 2 lists 
the upper limits we have used for comparison with 
the X-ray afterglows. This list has been compiled 
by using the strongest upper limits available for 
a specific burst. Also, for the more recent bursts, 
currently only reports in GCN and IAU Circulars 
exist, which should be considered temporary (this 
of course may be true for several of the older bursts 
as well). 

2.3. Predicting the optical flux by extrap- 
olating the X-ray afterglows 

We extrapolated the X-ray afterglow to the 
epoch and wavelength of the available optical and 
near-IR upper limits, as follows: we first obtained 
a range of power-law indices p of the electron en- 
ergy distribution from the la extremes of the X- 
ray power law temporal and spectral indices. To 
do this, we assumed eight cases in the standard 
afterglow model, which were derived from a com- 
bination of the following: the cooling frequency 
v c was above or below the X-rays, the epoch of 
the X-ray observation was before or after the jet- 
break in the light curve, and the fireball was ex- 
panding into a constant or a stellar wind density 
medium. For convenience, we have listed these 
combinations and the dependence of 8 and (3 on 
p in Table 3. Some of these cases were ruled out 
by the fact that the two electron, inferred from 
the X-ray 8 and (3 values respectively, are incom- 
patible with each other. We did not, however, 
limit the allowed range of p to p > 2 (for p < 2, 
one needs to put a high-energy cut-off on the elec- 
tron distribution; see e.g. Panaitescu & Kumar 
2001 or Bhattacharya 2001). From the obtained 
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values for p, we extrapolated in each case the 
+lcr and —1(7 X-ray flux to the epoch and fre- 
quency of the optical/near-IR upper limit. For 
this extrapolation, we chose the extreme possibili- 
ties, such as placing the cooling frequency at either 
the optical frequency or at the X-ray frequency, 
or the epoch of the jet-break right before or af- 
ter the epoch of the X-ray detection. In this way, 
we obtain the maximum and minimum values at 
optical/near-IR frequencies of the extrapolated X- 
ray afterglow. Sometimes such an extrapolation is 
ruled out by the standard afterglow model, and 
the result was ignored. This happens for exam- 
ple when v c reaches when extrapolating back 
in time, while p was derived from the case where 
Vy. < v c . 

By choosing the lowest and highest extra- 
polated flux, we obtain the flux range in which 
one would expect an optical counterpart to be 
found. Any upper limit above this range cannot 
be used to constrain the brightness of the optical 
afterglow. Any upper limit below this flux range 
will define this burst as physically 'dark'. Upper 
limits found within the flux range will be called 
'grey', and are subject to further study. Of course, 
the optical afterglow of grey bursts could very well 
be below the extrapolation of the X-ray afterglow, 
but so far there is no observational evidence sup- 
porting this assumption. 

In some cases, the electron index p inferred from 
the available spectral and temporal index reached 
a value below 2. Since this index is generally found 
to be > 2 for detected optical afterglows (but see 
some of the results of Panaitescu & Kumar (2001), 
who find indices of p w 1.5 for some afterglows), 
we have also calculated the lowest possible extra- 
polation with p fixed at 2. In most cases, however, 
the p = 2 extrapolation is well below the lowest 
minimum extrapolation from the available indices. 

In the case of GRB 040223, a potential dark 
burst, the temporal and spectral indices of the 
X-ray afterglow (Ticngo et al. 2004) are incom- 
patible with each other within the fireball model 
used here. We have listed the X-ray properties in 
Table 1, but have not extrapolated the X-ray flux 
to the optical waveband. 



3. Results 

In Figure I we show the results of the X-ray ex- 
trapolations and the various available upper limits 
at optical and near-IR wavelengths. For each X- 
ray afterglow, we only show the strongest optical 
and/or near-IR upper limit, as well as the mini- 
mum and maximum temporal extrapolation of the 

— 1(7 and +1(7 X-ray flux to the frequency of the 
upper limit (solid lines). The dashed line is the 
extrapolation in the case p = 2, obtained for the 

— 1(7 X-ray flux and the minimum temporal and 
spectral extrapolation. 

There are only 3 GRBs out of the sample of 20 
that we have found to be physically dark in this 
classification scheme: GRB 970828, 000210 and 
001025A. Two of these bursts have already been 
noted to be dark by extrapolation of the available 
X-ray measurements (GRB 970828, Groot et al. 
1998, and GRB 000210, Piro et al. 2002). 

In Table 4, we have indicated the dark bursts, 
together with the minimum amount of extinction 
(expressed in magnitudes) needed in the observer 
frame to explain the non-detection of the source. 

We have also plotted i?-band magnitudes for 
all available afterglows together with all i?-band 
upper limits in a magnitude-versus-time diagram 
(time being the delay between the GRB trigger 
and the observation; Figure 2), and indicated 
which upper limits belong to the dark bursts. Up- 
per limits in other bands with an available X-ray 
spectral index were converted to i?-band magni- 
tude following the same method as used in the X- 
ray extrapolations (i.e., the resulting i?-band up- 
per limit is the extrapolation with the least con- 
straining magnitude). We have drawn a power law 
with temporal index a = —1.5 passing through 
25th magnitude at 2 days, which separates the 
three dark bursts from practically all detected af- 
terglows. This power law is not based on any 
particular physical model, but merely serves as a 
guide in the diagram. 

The validity of this division, particularly at 
early times, has to be investigated, but at the cur- 
rent moment the lack of early X-ray afterglows 
prevents the study of dark bursts at early times 
in two ways: 1) early X-ray afterglows will gener- 
ally provide more accurate positions than the ini- 
tial GRB position, which gives larger telescopes, 
which often have a smaller field-of-view, the pos- 
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sibility to concentrate on deep searches, and 2) 
early X-ray afterglows will prevent the increase in 
the error when extrapolating the X-ray afterglow 
in the time-domain to its optical/near-IR upper 
limit. This lack will be amended with the early 
afterglows available from Swift. 

Two bursts with an optical afterglow fall be- 
low this line, GRB 020322 and GRB 030115. For 
GRB 020322, the afterglow was detected at R = 
23.26 (Bloom et al. 2002b) after the identification 
of the X-ray afterglow with XMM-Newton (Ehle 
et al. 2002). We have extrapolated the X-ray af- 
terglow for this burst in the same way to its ex- 
pected optical detection (see Figure lz), and we 
find it to be compatible within the la range of the 
X-ray extrapolation. GRB 030115 was R = 24.3 
one day after the burst and is very red, with an 
R — K color of about 6 mag (Levan et. al., ApJ 
submitted). Therefore, it is likely that the burst 
suffers from large extinction. Unfortunately, no 
X-ray afterglow has been reported for this burst. 

The other upper limits below this line belong 
to (in time order) GRB 030324 (Moran et al. 
2003) (first two upper limits), GRB 021113 (Levan 
et al. 2002), GRB 030324 (Moran et al. 2003), 
GRB 000615 (Stanek et al. 2000), GRB 990704 
(Rol et al. 1999), GRB 000214 (this paper), 
GRB 001109 (Vreeswijk et al. 2000), GRB 981226 
(Lindgren et al. 1999), GRB 990806 (Greiner et al. 
1999, equivalent i?-band magnitude), GRB 991105 
(Palazzi et al. 1999), GRB 981226 (Bloom et al. 
1998, equivalent i?-band magnitude) and GRB 000830 
(Jensen et al. 2000). Seven of these bursts 
(GRB 981226, GRB 990704, GRB 990806, GRB 991105, 
GRB 000214, GRB 001109 and GRB 000615) are 
included in our sample of analyzed afterglows, and 
shown in Figure 1. For two of these (GRB 000214 
and GRB 001109), the optical/near-IR limit is 
very close to the minimum extrapolation from 
X-rays, putting them at the border of physically 
defined grey/dark bursts. In fact, as this mini- 
mum extrapolation is an extreme possibility for 
the afterglow behavior, it would not be unreason- 
able to qualify GRB 000214 and GRB 001109 as 
dark bursts. However, to adhere to our definition 
of physically dark bursts, we choose not do so 
here. 



4. Discussion 

Figure 2 shows that in an i?-magnitude-versus- 
time diagram, the physicallly defined dark bursts 
are located in a distinct region at the faint end 
of optically detected bursts. The detection of 
GRB 020322 and GRB 030115, whose magnitudes 
are comparable to the upper limits of the two dark 
bursts GRB 000210 and GRB001025A, suggests 
that dark and detected bursts form a continuous 
group, rather than two clearly separated samples, 
and that the afterglows of most of the dark burst 
candidates might be found by deeper searches. 
Such afterglows may very well be highly reddened, 
as is the case for GRB 030115, and would require 
searches to be performed in the near-IR rather 
than in the optical wavebands. One would ex- 
pect that some dark bursts occupy the region of 
detected bursts; unfortunately, the X-ray data do 
not allow us to constrain the upper limits in this 
sample enough to find dark bursts in this region. 
Specifically, a well defined spectral index will con- 
strain the result of the X-ray extrapolation to 
optical/near-IR wavelengths. Such a constrained 
X-ray extrapolation, combined with an optical up- 
per limit, can show if physically dark bursts also 
occur in the upper part of Figure 2. More signif- 
icantly, our method can be used as a diagnostic 
tool for the overall determination of dark bursts, 
provided that early X-ray and optical observations 
are available, to compare the calculation of the 
extrapolated values to the observed limits. 

Here, we speculate on three possible causes for 
the (physical) darkness: high-rcdshift bursts, dif- 
ferent afterglow properties and extinction within 
the host galaxy. However, we wish to note that 
our sample of three bursts is too small to draw 
any definite conclusions. 

High-redshift bursts will go undetected at op- 
tical wavelengths due to extinction by the Lya 
forest, although this would require an X-ray af- 
terglow which is intrinsically much harder, as its 
redshifted flux is still detected in the 2-10 keV 
X-ray range. Furthermore, two dark bursts have 
a redshift from the detection of their host galaxy 
(GRB 970828, z = 0.9578 (Djorgovski et al. 2001), 
and GRB 000210, z = 0.846 (Piro et al. 2002)), 
which is well within the sample of detected after- 
glows and does not suggest a high-rcdshift ori- 
gin for dark bursts. The redshifts of the two 
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detected bursts with very dim afterglows, which 
fall within our region of dark bursts (GRB 020322 
and GRB 030115), are not known. However, for 
GRB 030115 the possible host is visible at R = 
24.5 (Garnavich 2003), putting its redshift at z < 
4.5. 

Dark bursts could form a subset of GRB af- 
terglows with different afterglow properties. De 
Pasquale et al. (2003) have found that the X-ray 
fluxes of dark bursts are different than the X-ray 
fluxes of optically detected bursts. Their sam- 
ple of dark bursts, however, is not physically de- 
fined as in this current study, but uses the obser- 
vational definition, which includes non-detections 
that might be attributed to observing conditions. 

To see if there is anything remarkable in the 
prompt gamma-ray properties, we have compared 
the 40-700 keV gamma-ray fluences of the three 
phsyically dark bursts with a general set of bursts, 
obtained from the BeppoSAX GRBM data and 
the BATSE data (the BATSE 40-700 keV fluences 
were obtained from spectral fits to the BATSE 
archival data, while the GRBM data were ob- 
tained from literature). This comparison set was 
selected by choosing bursts for which narrow error 
boxes are available and have had optical follow-up 
campaigns, to make the comparison set similar to 
the set of GRB data of the physically dark bursts. 
The set includes among others all the grey bursts. 
For GRB001025A, the Ulysses 25-100 keV fluence 
was converted to the 40-700 keV range, using a 
conversion factor obtained by converting a set of 
28 BATSE burst templates with known spectral 
shapes (cf. Bloom et al. 2001). The result is shown 
at the top of Figure 3. No obvious trend is visi- 
ble, although we note that the dark bursts appear 
to have a higher than average 40-700 keV fluence. 
This is also noticeable in the bottom figure, where 
we have plotted the peak fluxes of the dark bursts 
together with those for a general sample of bursts: 
the three dark bursts are in the higher part of the 
sample. However, our set is too small to meaning- 
fully apply any formal statistical test, and we note 
it here as a possible trend. 

Such a correlation between the peak flux and 
the faintness of the optical afterglow could happen 
in the case of an external scenario for the GRB and 
its afterglow: if the prompt gamma-ray emission 
is the result of the interaction with an external 
medium instead of internal shocks, a higher den- 



sity of the medium would then result in a higher 
gamma-ray peak flux of the burst. A higher den- 
sity would, however, also result in a higher ab- 
sorption depth (provided the path length is long 
enough) and the optical/UV emission would go 
undetected. 

Extinction as a likely cause for dark bursts 
was already suggested by Groot et al. (1998) for 
GRB 970828. More recently, Klose et al. (2003) 
explored the available AT-band upper limits of 
GRBs and also found evidence for extinction to 
explain dark bursts. The two dark bursts with 
an identified host galaxy have already been dis- 
cussed before. GRB 970828 has one of the deep- 
est early upper limits available, as well as a host 
galaxy and a redshift. Djorgovski et al. (2001) 
show that a possible dust lane or giant molecu- 
lar cloud could provide the amount of extinction 
needed to obscure the afterglow. For GRB 000210, 
the X-ray spectrum shows significant Ah absorp- 
tion; the amount of dust needed to obscure an op- 
tical afterglow at the redshift of the host galaxy is 
consistent with this column density. Both findings 
favor dust extinction in the host galaxy for these 
two dark bursts. 

Extinction should also be noticeable in opti- 
cally detected bursts, and the R-K color of 6 
mag for GRB 030115 indeed suggests that this 
burst is heavily extincted. However, only very few 
bursts indeed show large extinction in the host 
galaxy. Even some notably faint afterglows, like 
GRB 980613 (Hjorth et al. 2002) and GRB 000630 
(Fynbo et al. 2001), have a low (< 0.5 magni- 
tude) extinction value in the i?-band (in addition 
to galactic extinction). It is therefore quite pos- 
sible that we are simply not finding many heav- 
ily extincted bursts because these are still below 
the generally reported detection limits. However, 
the last 10 optical/near-IR follow-up campaigns of 
HETE-2 SXC error boxes resulted in 8 detected 
afterglows. This suggests that for the majority 
of the 60% of the cases where no afterglow was 
found for a BeppoSAX GRB, this was due to ad- 
verse observing conditions. Our sample and the 
SXC sample suggest a fraction of physically dark 
bursts more like 10-20%. This may indicate that 
the few dark bursts we find here, are indicative of 
the size of the actual (physical) dark burst sample. 
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5. Conclusions 

The extrapolation from the X-ray afterglow to 
optical/near-IR wavelengths shows that only very 
few bursts deserve the designation 'dark burst'. 
The optical upper limits obtained for these dark 
bursts are within a region at the fainter end of 
the optically detected bursts, with two optical af- 
terglows found within this region. One of these 
two optical afterglows is highly reddened, which 
suggests that dark bursts can indeed be detected, 
albeit at very faint magnitudes. Detection of such 
dark bursts provides the necessary insights into 
the cause of their darkness. 

Our method of extrapolation provides a use- 
ful diagnostic tool for setting a lower limit on the 
expected magnitude of the optical counterpart. 
When an early X-ray light curve and spectrum 
are available, for example from the XRT on board 
Swift, one can set a 'dark burst upper limit'. Any 
afterglow not detected above this upper limit is 
a potential candidate for a high rcdshift or an ex- 
tinguished afterglow, or possibly even an afterglow 
produced by different physical mechanisms. Set- 
ting a search area using the currently available af- 
terglow upper limits, as in Figure 2, gives future 
afterglow searches a tool for picking out such po- 
tential dark bursts even without an X-ray after- 
glow. 

It is, however, not possible to draw definite con- 
clusions on the cause of darkness, due to the small 
set of dark bursts we obtained. We note that in the 
gamma-ray properties, dark bursts appear to have 
higher than average fluences and peak fluxes, but 
we cannot statistically quantify this result yet. We 
also find that at least one burst showing large red- 
dening (GRB 030115), is within the same region of 
magnitude-versus-time diagram (Figure 2) as the 
dark bursts. 

Deep upper limits within one day, as well as 
early localizations and rapid follow-up observa- 
tions will complete and enhance Figure 2. The 
Swift mission can certainly help filling the gap 
for early times: it will provide early X-ray after- 
glow spectra and light curves, as well as rapid and 
precise localizations, needed for early (robotic) 
searches and deep searches with larger telescopes. 
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Fig. 1. — The ranges obtained for the optical counterpart fluxes for all GRBs in Table 1, by extrapolation of 
their X-ray measurements to the epoch and frequency of their optical upper limits; these limits are indicated 
by a downward-pointing arrow. The solid lines show the two extreme extrapolations, within which an optical 
counterpart should be found: the bottom line was obtained by extrapolating the — la X-ray flux according to 
the lowest possible temporal and spectral slopes, and the top line was obtained by extrapolating the + la X- 
ray flux using the highest possible temporal and spectral slopes. The dashed lines follow from extrapolating 
— la X-ray flux with the lowest temporal and spectral slope with p = 2. 
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Fig. 2. — An R-magnitude-versus-time diagram, containing all the i?-band upper limits (asterisks) from 
Table 2, as well as all known detected afterglows (diamonds). For most afterglows, we have taken the first 
available i?-band data point. The dark bursts listed in table 4 are indicated by the three filled circles. 
The power law with a = —1.5 indicates the region that separates the dark bursts from almost all detected 
afterglows. 
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Table 1 

Overview of X-ray afterglows of bursts with optical non-detections. We have also 
indicated the minimum and maximum values of the electron index p, which correspond to 
the minimum and maximum extrapolations in figure 1 respectively. 
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a Optically detected burst. See Sect. 3 for details. 

b The values for 5 and (3 are incompatible with each other within the variants of the fireball model used in 
this paper, yielding no possible range for p 
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Table 2 

Overview of optical/near-IR upper limits with X-ray afterglow data available 
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20.90 


0.04 


Greiner et al. (2000) 


GRB 001109 


K 


0.400 


19.90 


0.04 


Castro Ceron et al. (2004) 


GRB 010214 


R 


0.893 


23.50 


0.05 


Rol et al. (2001) 


GRB 010220 


R 


0.361 


23.50 


0.85 


Berger et al. (2001) 


GRB 010220 


K 


0.974 


17.10 


0.85 


Li Causi et al. (2001) 


GRB 011030 


R 


0.300 


21.00 


0.39 


Mohan et al. (2001) 


GRB 011030 


K 


0.520 


19.00 


0.39 


Klosc et al. (2001) 


GRB 011030 


R 


1.000 


23.61 


0.39 


Rhoads et al. (2001) 


GRB 011030 


V 


42.905 


27.6 


0.39 


Levan et al. (2005) 


GRB 020321 


R 


1.099 


23.60 


0.10 


Salamanca et al. (2002) 


GRB 020322 


R 


0.329 


23.26 


0.08 


Bloom ct al. (2002b) a 


GRB 020427 


V 


44.603 


28.7 


0.03 


Levan et al. (2005) 



a Optically detected burst. See Sect. 3 for details. 
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Table 3 

The temporal (6) and spectral (f3) indices as a function of the electron power-law 

DISTRIBUTION INDEX p IN THE SLOW COOLING CASE, FOR EIGHT DIFFERENT SCENARIOS. SEE E.G. SARI 

ET AL. (1998) and Chevalier & Li (1999). 



Sip) 


Pip) 


Sip) Pip) 


v > v c 




v < v c 


Pre jet-break (isotropic) 






ISM (p = constant) (3p - 2) /4 
Stellar wind (p cx R~ 2 ) (3p - 2)/4 


p/2 
p/2 


3(p-l)/4 (p-l)/2 
(3p-l)/4 (p-l)/2 


Post jet-break (collimated) 






ISM (p — constant) p 
Stellar wind (p oc R~ 2 ) p 


p/2 
p/2 


p (p - l)/2 
p (p - l)/2 



Note. — We have chosen (5 and /3 in such a way that they are positive in optical and X-rays for positive 
p, i.e. the afterglow flux follows F(i>,t) oc v~Pt~ 5 . 



Table 4 

Burst classification and amount of (host galaxy) extinction inferred (in the observer 

frame) 



burst 


filters 


classification 


extinction (mag) 


GRB970402 


R 


grey 




GRB970828 


R 


dark 


> 4.304 


GRB971227 


R 


grey 




GRB981226 


R 


grey 




GRB990704 


R 


grey 




GRB990806 


B 


grey 




GRB991014 


K, R 


grey 




GRB000210 


R 


dark 


> 1.708 


GRB000214 


K, R 


grey 




GRB000528 


R 


grey 




GRB000615 


H, R 


grey 




GRB001025A 


R 


dark 


> 0.690 


GRB001109 


K, R 


grey 




GRB010214 


R 


grey 




GRB010220 


K, R 


grey 




GRB020321 


K, R 


grey 
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